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Abstract: Moisture equilibrium data of Urtica urens leaves by desorption and adsorption were 
determined at 30, 40, and 50°C. The isotherms were determined by a static gravimetric method for 
various temperature and humidity conditions. Five mathematical models were used to fit 
the experimental data. A nonlinear least-squares regression program was used to evaluate the constants 
of the five desorption and adsorption isotherm models GAB and Modified Halsey models were adequate 
to describe the sorption characteristics of the samples. Isosterie heats of desorption and adsorption were 
calculated by applying the Clausius-Clapeyron equation to the sorption isotherms at different 
temperatures. It decreased with increasing moisture content. We recognized the linear relation exists 
between the enthalpy and entropy of the sorption reaction. 
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I. INTRODUCTION 

The homeopathic remedy Urtica urens leaves is prepared using the entire herb known as annual 
nettle, while it is in blossom and is used to treat a number of health conditions, especially nettle rash or hives. 
Commonly known as dwarf nettle, annual nettle or small nettle, Urtica urens is a semi-woody plant belonging 
to genus Urtica that grows annually. Traditionally, people have been using this herb as a medication for 
skincare as well as alleviate skin inflammation or itchiness. In effect, scientists have especially studied the herb 
Urtica urens for its potential in supporting the skin when external irritants result in certain complaints [1]. 

The medicinal and aromatic plants have a great importance for both the pharmaceutical industry and 
the traditional. The stability of a dehydrated medicinal plant is influenced by its water activity. This stability is 
mainly a consequence of the relationship between the equilibrium moisture content and its corresponding water 
activity. Knowledge of the sorption properties of foods is of great importance in medicinal plant, especially in 
the quantitative approach to the prediction of the shelf life of dried medicinal plant. Equations for modeling 
water sorption isotherms are of special interest for many aspects of medicinal plant preservation by 
dehydration. 

In this study, temperatures of 30°C, 40°C, and 50°C were chosen, to typify tropical storage 
conditions. Using an experimental approach, the equilibrium curves are determined by the saturated salt 
solution method. The experimental sorption curves are described by five different models to identify the most 
appropriate mathematical model for a better description of the product equilibrium state: GAB's model, 
modified Henderson's model, modified Halsey's model, modified Oswin's model. Then, we investigate from 
the experimental data the water activity optimal for the storage of Urtica urens. 

Further analysis of sorption isotherm data by application of thermodynamic principals can provide 
information regarding differential enthalpy or the isosterie heat of sorption and differential entropy. 
Differential enthalpy of sorption gives an indication of the amount of bound water existing in the product. The 
differential entropy of material is proportional to the number of its available sorption sites at a specific energy 
level [2]. Subsequently we determine the thermodynamic functions (differential enthalpy, differential entropy). 
Moreover, we recognized the linear relation exists between the enthalpy and entropy of the sorption reaction. 

II. Materials And Methods 

2.1 Experimental procedure 

The plant Urtica urens used in our study was collected at Tahanouat, Marrakesh Morocco. The 
hygroscopic equilibrium could be achieved by a dynamic or static method. In the present work, a static method 
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is used [3]. The method is based on the use of saturated salt solutions to maintain a fixed relative humidity Rh. 
The mass transfers between the product and the ambient air are assured by natural diffusion of the water 
vapour. The atmosphere surrounding the product has a fixed relative humidity for every working temperature 
imposed on the system. 

Six salts were chosen { KOH, (MgCl 2 , 6H 2 0), K 2 C0 3 , NaN0 3 , KC1 and (BaCl 2 , 2H 2 0)} so as to have 
a range of water activity of 0.0572-0.0898 [4]. The experimental apparatus consisted of six glass jars of 1 liter each with 
an insulted lid. Every glass jar was filled to quarter depth with a saturated salt solution. Duplicated samples each of 0.4 g 
(+0.0001g) for desorption and 0.1 g (+0.0001g) for adsorption were weighed and placed into the glass jars. The 
weight recording period was about 3 days. This procedure continued until the weight was constant. The 
equilibrium moisture content of each sample was determined in a drying oven at 105 °C for 24 h. The 
hygroscopic equilibrium Urtica urens leaves was reached in then days for desorption and eighth days for 
adsorption. 



2.2 Modelling equations 

The relationship between equilibrium moisture content, water activity, and temperature for Urtica 
urens leaves was predicted in literature by several mathematical models [5]. In our study five models have been 
applied: GAB's model, modified Henderson's model, modified Chung-Pfost model, modified Halsey model, 
and modified Oswin's model. The equations expressing those models are given below in table 1. X eq represents 
equilibrium moisture content (%MS); a w water activity; A, B, and C coefficients, dependant to the temperature 
or not, to determine by smoothing the experimental curves of sorption and 0 the temperature in °C. The 
correlation coefficient (r) was one of the primary criteria for selecting the best equation to fit the four models to 
the experimental data. In addition to r, the statistical parameters mean relative error (MRE) as a percentage 
was used to determine the quality of the fit [6]. Levenberg-Marquardt nonlinear optimization method using 
appropriate software is used for the calculation of model coefficients that describe the equilibrium curves and 
their statistical parameters: the correlation coefficient (r) and mean relative error (MRE). 



X(Xeq ipred -Xeq iexp )- 

X( Xe£ li.exp-X^q,exp) 2 



lOO N Xeq 
MRE = X 



(1) 



Where: 

Xeq. 

Xc q . pred 

N 
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i* Experimental 
i lh Predicted moisture 



number of data points, 
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Table 1: Mathematical models used in our study 



Expression of the model 



[1-Ba„] [l-Ba w +BCa w ] 



Modified Henderson 
Modified Chung-Pfost 
Modified Halsey 

Modified Oswin 



-ln(l-a w ) 
A(6> + B) 



- exp(-CXeq) 



x = (-exp(A + B.6Q 
eq [ ln(a w ) 

X eo =(A + B.6>) 



[9] 
[10] 
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2.3 Determination of the optimum conditions for storage 

The study of the sorption isotherms enables us to know the optimal relative equilibrium moisture for 
the conservation of a product as well as the water content of equilibrium to reach at the end of the drying. Also, 
it affords users with accurate information on how to handle a product during storage and conservation [12]. 
For this purpose, the optimal water activities of conservation (a wop ) were dogged. The whole of the 
experimental points is gathered on the same graph. Then, we model the isotherm of sorption by a polynomial 
equation of the third degree. The part closes to horizon tality corresponds to the zone of better stability of the 
product. This process of calculation consists in making a polynomial decomposition of the equilibrium 
moisture content X eq , for the whole of the experimental results, according to the water activity. This makes it 
possible to calculate the value for which the derivative second of X eq is cancelled "inflection point" and 
consequently optimal relative humidity for conservation. 

2.4 Differential enthalpy and differential entropy of sorption 

The differential enthalpy or isosteric heat of sorption ( AH d ) indicates the state of absorbed water by 
the solid material. The net isosteric heat of sorption ( Ah d ) represents the quantity of energy exceeding the 
heat of vaporization of water ( AH ) associated with the sorption process. 

Ah d = AH d - AH vap (3) 
The net isosteric heat of sorption can be calculated from the experimental data using the Clausius- 
Clapeyron equation [13]: 

r ^QnQ ] = -Ah^ (4) 
[d(l/T)\ x ^ R 

This relation requires determining the isotherms at various temperatures in order to calculate the 
variation logarithmic curve of the water activity according to the reverse of the temperature, for a water content 
fixed. However, starting from the modified equations, it is possible to obtain the analytical expression of the 
heat of sorption Ah d . This approach assumes that isosteric heat AH d does not vary with the temperature 
[14]. The relationship between the net isosteric heat Ah d and the differential entropy AS d of sorption is given 
by: 

Ah, AS ri 

(-lna w ) x = — d - & - (5) 

w x " R.T R 

The experimental sorption isotherm data was plotting in the form of ln(a w ) versus J_ for specific 

T 

moisture content (X eq ). We determined Ah d from the slope ( ^Kl \ and AS d from the intercept! ^ d \ . 

I R J I R J 

Applying this at different moisture contents allowed determining the dependence on Ah d and AS d . 

2.5 Enthalpy-Entropy Compensation Theory 

According to the theory of compensation, the linear relationship between the enthalpy Ah d and the 
entropy AS d for a specific reaction are given by: 

Ah d = Tp.AS d + AGp (6) 
The isokinetic temperature ( T ^ ) is a characteristic property of the food surface whose dimension is 
absolute temperature. It represents the temperature to which all the reactions of the series proceed to the same 
rate. The free energy ( AG ^ ) provides a criterion to evaluate if the process of the water sorption is 
spontaneous ( —AGp ) or not ( +AG^ ). This process is repeated for several values of equilibrium relative 
humidity with an aim of underlining the relation between the differential entropy of sorption and the 
equilibrium moisture content of the studied product. The isokinetic temperature ( 7^ ) and constant ( AGp ) 
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were calculated using linear regression. The isokinetic temperature Tp represents the slope of the linear 
function between the differential entropy and the differential enthalpy of sorption [2]. 

III. Resultas And Discussion 

3.1 Adsorption and desorption isotherms 

The equilibrium moisture content increases with decreasing temperature at constant water activity 
(Fig 1). The sorption isotherms present the characteristic S-shaped curve (Type II), typical of sorption 
isotherms of many plants and food materials [3, 4, 5]. 
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Fig. 1: Sorption isotherms of Urtica u 
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The presentation of adsorption and desorption experimental data in the same plot show that the plant 
does not have the same hygroscopic equilibrium behaviour (fig 2). The phenomenon of desorption-adsorption 
is irreversible. For the same constant relative moisture, the water content of desorption is higher than that of 
adsorption. Indeed, there exists a phenomenon of hysteresis. This characteristic is visible for Urtica urens 
plant. 

70 -i . 




Water actvitya w (-) 
Fig.2: Sorption hysteresis phenomenon of Urtica urens leaves 



Several hypotheses have been put to explain hysteresis. One of it is an analogy with a sponge. When 
the pores do not contain any more water, adsorption does not make it possible to recover the totality of initial 
water since this one imprisons air in the pores. Another hypothesis states that hysteresis in the sorption 
isotherm is a consequence of variation in the fraction of bound water present in the adsorption and desorption 
processes. The bound fraction being always larger on desorption than on adsorption [15]. 

3.2 Modelling of sorption experimental data 

The results of non linear regression analysis of fitting the sorption equations to experimental data of 
Urtica urens leaves was at three temperatures are presented in Tables 2 and 3. For all the models tested, 
parameters A, B and C are found to be temperature dependent. For desorption isotherms, GAB model is found 
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to be the best estimator for predicting the equilibrium moisture of the Urtica urens leaves. This model gives the 
highest correlation coefficient (r) of 0.9989 and mean relative error MRE (2.7865 %) at 8 = 40 °C. For 
adsorption isotherms, modified Halsey model is found to be the best estimator for predicting the equilibrium 
moisture of Urtica urens leaves. This models gives the highest (r) of 0.9918 and the MRE (2.8321 %) at 
9 = 50°C. 



Table 2: Estimated model coefficients, r, MRE of five equation model fitted to desorption isotherm of Urtica 
urens leaves 



Model's name 


9 (°C) 


A 


B 


c 


r 


MRE (%) 


GAB 


30°C 


0.9365 


151.7891 


5.5773 


0.9986 


2.1527 


40°C 


0.9803 


46.7567 


7.3417 


0.9989 


2.7865 


50°C 


0.9485 


112.8598 


5.6501 


0.9975 


3.0121 


Modified Henderson 


30°C 


0.0039 


-17.4841 


1.1021 


0.9697 


4.3719 


40°C 


0.0023 


-15.9677 


1.0599 


0.9791 


6.7184 


50°C 


0.0018 


-19.4031 


1.0534 


0.9775 


5.6338 


Modified Chung-Pfost 


30°C 


124.4991 


8.49 


0.1105 


0.9472 


8.7143 


40°C 


124.0871 


9.6501 


0.0633 


0.9411 


15.2637 


50°C 


146.5175 


-3.1678 


0.1068 


0.9492 


13.6218 


Modified Halsey 


30°C 


2.6725 


0.0274 


1.6101 


0.9955 


3.3821 


40°C 


2.3441 


0.0241 


1.5388 


0.9966 


3.2581 


50°C 


1.1786 


0.0206 


1.4927 


0.9965 


4.8321 


Oswin Modifie 


30°C 


4.2552 


0.2352 


1.9592 


0.9872 


5.3127 


40°C 


1.0984 


0.2525 


1.8779 


0.9915 


4.1287 


50°C 


0.2431 


0.2211 


1.8358 


0.9905 


7.7385 



Table 3: Estimated model coefficients, r, MRE of five equation model fitted to adsorption isotherm of Urtica 
urens leaves 



Model's name 


9(°C) 


A 


B 


C 




MRE( %) 


GAB 


30°C 


0.9538 


756.4351 


4.3744 


0.9912 


3.1829 


40°C 


0.9351 


134.8204 


4.7111 


0.9866 


3.1815 


50°C 


0.9237 


282.8957 


4.7116 


0.9915 


2.1923 


Modified Henderson 


30°C 


0.0034 


-10.9261 


1.0657 


0.9520 


6.8457 


40°C 


0.0019 


-12.1797 


1.1341 


0.9594 


8.1745 


50°C 


0.0006 


-9.0857 


1.2814 


0.9605 


5.7332 


Modified Chung-Pfost 


30°C 


112.5638 


5.186 


0.13 


0.9302 


13.7184 


40°C 


136.1907 


1.0462 


0.1357 


0.9402 


12.9637 


50°C 


153.7906 


-9.9718 


0.1523 


0.9481 


15.6013 


Modified Halsey 


30°C 


2.2536 


0.0193 


1.5443 


0.9881 


3.0821 


40°C 


2.4539 


0.0178 


1.5964 


0.9847 


3.3586 


50°C 


2.46 


0.0182 


1.6861 


0.9918 


2.8321 


Modified Oswin 


30°C 


1.9962 


0.2422 


1.8809 


0.9769 


4.3167 


40°C 


0.8772 


0.2146 


1.9685 


0.9774 


6.6218 


50°C 


0.7391 


0.177 


2.1286 


0.9819 


8.6253 



3.3 Measurement of optimal water activity of Urtica urens leaves 

A quality product deserves quality treatment. The heat and process engineering specialists must 
provide users with accurate information on how to handle a product during storage and conservation. For this 
purpose, the optimal water activities of conservation ( a w ) were determined. The sorption isotherm curve can 
be described as a polynomial equation of the third degree, the central part or "stage" is the best area of product 
stability (fig 3). 
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'3* 0,0 0,2 0,4 0,6 0 : 8 1,0 
rt Water activity a (-) 



Fig.3: Determination of optimal water activity for conservation of Urtica urens leaves 

This calculation method consists of decomposing polynomial of equilibrium moisture content Xeq, 
for all the experimental results for each product based on water activity (equation 7). The values of the optimal 
water activity conservation of Urtica urens leaves (fl = 0.3532) are in perfect agreement with the results 
obtained for other products [16, 17]. The optimal water activities are ranged in (0.3-0.4). 
X eq =1.6442+57.7527a w -167.6854a w 2 + 157.38322a w 3 (7) 

3.4 Enthalpy and entropy of sorption 

The net isosteric heat of sorption (Ah d ) values are calculated from the equilibrium moisture data at 
different temperatures using Clausius-clapeyron equation (4). The variation of the heats sorption of Urtica 
urens leaves with equilibrium moisture content is shown in figure 4. 

The net isosteric heat of desorption an adsorption of Urtica urens leaves can be expressed 
mathematically as a polynomial function of equilibrium moisture content: 
Desorption : Ah s =l 1.753 - 1.2126 X eq + 0.0436 X e 2 q - 0.0005 X^ q (8) 
Adsorption : Ah s = 19.8186 - 3.5446 X eq + 0.21 16M e 2 q +0.004 X^ q (9) 

The differential entropy (AS d ) is plotted as a function of moisture content in figure 4. Once again, 
the entropy data display a strong dependence on moisture content. The experimental desorption and adsorption 
data, respectively, conform to polynomial relation as represented by the equations: 
Desorption : AS d = 40.0995 - 4.2976 X eq + 0. 1727 X e 2 q - 0.0023 X^ q (10) 

Adsorption : AS d = 62.5820 - 10.9222 X eq + 0.6521 X 2 q - 0.0125 X^ q (11) 




Equilibrium moisture contentXeq (.b) Equilibnumn moisture content X ^ (.b) 



Fig.4: Net isosteric heat and differential entropy of sorption of Urtica urens leaves as a function of equilibrium 
moisture content 
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3.5 Enthalpy-entropy compensation theory 

The Ah d and AS d values for adsorption and desorption, at given moisture contents, were calculated 
by linear regression equation using equation 5. It was assumed that, at specific moisture content, Ah d and 
AS d did not vary with temperature. Figure 5, shows a linear relation, with a correlation of determination (r 
=1). This indicates that the compensation exists. The parameters T p and AG p (Equation 6) were calculated 
from the data by linear regression, and the values are detailed in table 4 




4 6 8 10 12 



AS.(J.mol' : .K" : ) 
Fig.5: Compensation theory of enthalpy/entropy 



Table 4: Characteristic parameters for Ah d I AS d Relationship 





T„(K) 


AGpCLmol 1 ) 


desorption 


340.96 


-623.62 


adsorption 


308.08 


374.35 



IV. Conclusion 

It's now widely accepted that sorption data provide valuable information for industrial drying in order 
to preserve and store the seasonal plants and make them available to consumers all year round. The sorption 
data may be used not only for the identification of the optimum residual moisture content of the final product, 
but also for the calculation of the drying time of hygroscopic substances. 

The moisture sorption curves of Urtica urens leaves were experimentally investigated at three 
temperatures 30, 40 and 50 °C. The isotherms have a sigmoid shape (Type II) which is common for many 
hygroscopic products. The hysteresis phenomenon was distinctly observed. Among the sorption models chosen, 
the GAB model was the best model describing the equilibrium moisture data for desorption, and the modified 
Halsey model was the most suitable to estimate adsorption isotherms. The optimal water activities of 
conservation were determined. Net isosteric heat of sorption and differential entropy values of Urtica urens 
leaves were found to be a polynomial function of moisture content for desorption and adsorption. Enthalpy- 
entropy compensation theory could be successfully applied to the moisture sorption behaviour of Urtica urens 
leaves. 
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